Proposed Detection of Time Reversal Symmetry in Topological Surface States 
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By employing T-matrix approach, we investigate a nonmagnetic impurity located on the surface of three- 
dimensional topological insulators, where the time reversal symmetry is preserved. It is shown that the images 
of the local density of states (LDOS) around the single impurity have the dip-hump structures with six-fold 
symmetry at different bias voltages. The peaks are produced by quasiparticle interference while the local minima 
at the backscattering wave vectors are due to the absence of backscatterings. With increasing the bias voltage, 
the peaks and dips move forward to the location of the impurity. These dips at the backscattering wave vectors 
in the LDOS spectra can be regarded as a signature of the time reversal symmetry in the topological surface 
states, which could be observed by scanning tunneling microscopy. 

PACS numbers: 73.20.-r, 72.10.-d, 72.25.-b 



In recent years much attention has been focused on the 
topological surface states due to their potential applications in 
quantum computing or spintronics [1,2]. The Dirac-cone-like 
electronic states, existing on the surface of three-dimensional 
bulk insulating materials, have been observed in Bii_ x Sb x [3], 
Bi 2 Sb 3 [4], Sb 2 Te 3 [5], Bi 2 Te 3 [5,6], TlBiSe 2 andTlBiTe 2 [7], 
by angle-resolved photoemission spectroscopy (ARPES). The 
surface energy band structure was determined by employing 
k - p theory [8], where an unconventional hexagonal warping 
term plays a crucial role in fitting the ARPES data [3-7] and in 
explaining the scanning tunneling microscopy (STM) experi- 
ments [9-13]. 

It is known that the novel topological surface states possess 
time reversal symmetry (TRS), which leads to the absence of 
backscatterings. Such a symmetry induces unusual electron 
transport properties on the surface of topological insulators 
[9-13]. Therefore, it is interesting to explore the TRS in order 
to understand thoroughly the physical properties of topolog- 
ical insulators. In this work, we give a proposal for experi- 
mental detection of the TRS in the topological surface states. 
We note that the oscillations of local density of states (LDOS) 
induced by a step or line defect can be explained satisfically 
by quasiparticle interference (QPI) and seem not to exhibit 
the existence of the TRS due to no direct evidence showing 
the absence of backscatterings. [9-13]. Here we employ a 
nonmagnetic impurity located on the surface of topological 
insulators in order to detect this symmetry. Such a single im- 
purity has also been used to judge the superconducting order 
parameter symmetry in the cuprates and the FeAs-based su- 
perconductors [14-16]. In Refs. [17,18], the Fourier transfor- 
mation of the LDOS induced by a nonmagnetic impurity was 
calculated. However, the oscillations of the LDOS with the 
hexagonal warping term and the role of the TRS in real space 
are never discussed previously. 

The momentum space Hamiltonian describing the surface 
states of three-dimensional topological insulators reads [8] 

h o = Yj C k[(^~T ~ ^ + v ( k *Py ~ k y°~^ + ^k^JCk, (1) 



where = (c^, c^), / and Ci(i = x,y,z) are the 2 x 2 unit 
matrix and the Pauli matrices, respectively, m* is the effective 
mass of electrons, which is usually very large for the topolog- 
ical insulators, \i is the chemical potential to be determined 
by doping [10], v is the strength of the Rashba spin-orbit cou- 
pling, the last term is the so called hexagonal warping term, 
and 0k = k x (k 2 x - 3ky). Taking the linear transformations 
Ckr = Z 5=0,1 Mksaksfiks and c^i = 4= 2 s =o,i fik S aks+i ftks 
with a^ s = 1 - s + s(ik x + k y )k~ l , = s + (1 - s)(ik x - k y )k~ l , 
®ks = cos j- + (-l) s sirij-, 6^ = arctan^-, and k = Jk^ + k*, 
we have H = ^ks Eks^l^ks with the eigenenergies E ks = 

Now we investigate a singe impurity located at the ori- 
gin r = (0, 0) on the surface of three-dimensional topologi- 
cal insulators. Without loss of generality, such an impurity 
can be described by the Hamiltonian H{ mp = U^a C lo- CO(T + 
V m (<4^oT ~ c ii c ®0- ^ ere U and V m are the nonmagnetic and 
magnetic parts of the impurity strength, respectively. We solve 
the total Hamiltonian H - Hq + H[ mp by employing T-matrix 
approach [15,16] and finally obtain the analytical expression 
of the LDOS around the impurity 
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where N is the number of lattices, p^ s '(r,cj) represents the 
LDOS contributed by the scattering process (ks) — > (kV), 
and the Green's functions G^ s \ioj n ) = G® (ioj n )S^S SS ' + 
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G^ s (icL> n ) = (ico n - E ks )~ l an 



) with the bare Green's functions 
and the T matrix 
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FIG. 1: (Color online) The stereographs of the LDOS p(r, oj) around a unitary impurity, i.e. U — > oo, at different bias voltages. 
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Here, A(ioj n ) = ± when V m = 0, we have 

T~^ s (iaj n ) = 0. This means that there are no backscatter- 
ings at each constant-energy contour oj = E ks [17,18]. In 
other words, the surface states with opposite momentum and 
the same s are still orthogonal and do not mix after being 
scattered by a pure nonmagnetic impurity. So there is no 



interference between the two surface states connected by a 
backscattering wave vector. From Eq. (2), we can also see 
that the backscattering processes have no contribution to the 
LDOS, i.e. Pkf(r, oj) = 0. Therefore, the LDOS induced 
by a nonmagnetic impurity should has the local minima at 
the backscattering wave vectors. Such the dips are produced 
by the absence of backscatterings rather than the QPI, which 
could be used to detect the TRS in the topological surface 
states. 



According to Eq. (2), we can calculate the LDOS around 
the single impurity, which is measured by STM experiments. 
In our calculations, we choose a 120 x 120 A 2 lattice with N = 
1800 x 1800 and employ the physical parameters in Bi2Sb3, 
i.e. A = 250.0 eV-A 3 , v = 2.55 eV-A, and ji = 0.334 eV 
[8,10]. 

In Fig. 1, we present the stereographs of the LDOS induced 
by a unitary impurity (U — > oo) at oj = -120, 0, 120 (meV). 
Such the infinite potential has been realized by the Zn ions 
in the cuprates [14], which is also expected to exist on the 
surface of topological insulators. Obviously, p(r, oj) has a dip- 
hump structure with a six-fold symmetry at each bias voltage. 
When the bias voltage oj increases, these peaks and dips move 
forward to the location of the impurity. On the impurity site, 
we have p(r, 6L>)| r =(o,o) = 0. 

In order to understand clearly the energy-dependent peaks 
and dips in the LDOS, we depict the constant-energy con- 
tours of the surface energy band at oj = -120, 0, 120 (meV) 
in Fig. 2. There are six modulation wave vectors q^i, q^2, 
qKB, <1mi» Qm2, QMfi, and their corresponding symmetric wave 
vectors induced by the QPI. q KB and q M B are the backscatter- 
ing wave vectors connecting two opposite K points and two 
opposite M points, respectively. According to the symme- 
try of the constant-energy contours, we have |q#i| = ^Iq^sL 

lq;d = T^lq^l, IqMil = \\<{mbI and |q M 2l = %qMBl 
We also note that the wave vectors q/n, q^#, and qM2 (<1mi» 
qMB, and q^) are along the same modulation direction, i.e. 



qKil IqKBl IqMi <S1m\I I<1mbI I<\ki)- The an § le between the 
two modulation directions is integer • 60° + 30°. Due to the 
existence of a backscattering wave vector and two QPI wave 
vectors in the same direction, we shall observe below that the 
peaks and dips in the images of the LDOS seem not to locate 
exactly at these modulation wave vectors. When oj = -120 
meV, |q M2 | ~ l<tol- With increasing the energy, all the mod- 
ulation wave vectors become longer, but the increasing length 
of each wave vector is different. Therefore, the competition 
among the wave vectors produces abundant energy-dependent 
structures in the LDOS. The oscillation periods corresponding 
to these wave vectors have the relations: T K \ = = 2T KB , 

TK2 = \it\ = if K ^ Tm = \it\ = 2Tmb > and Tm = = 

jj=T M B- 

Fig. 3 shows the images of the LDOS with a moder- 
ate nonmagnetic potential U = 5.0 eV and V m = eV at 
oj = -120,0,120 (meV), where the TRS persists. When 
oj = -120 meV, the LDOS p(r,a>) has the local minima at 
r = T K2 ,T M 2,Tkb, and T MB . Because \T M2 \ « \T KB \ and 
\Tki\ ~ \TmbI these dips are produced by both the QPI and 
the absence of backscatterings. With increasing the bias volt- 
age, the local minima move forward to the origin and form a 
hexagonal sinkage around the impurity at oj = meV. Obvi- 
ously, p(r, oj) also has the peaks at r = T K \ and T M \ while the 
dips show up at 27^2 in Fig. 3(b). Because T K \ = 2T KB 
and T M \ = 2T M b, the effect of the QPI on the LDOS is 
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FIG. 2: (Color online) The constant-energy contours of the surface 
state band at different energies. The modulation wave vectors in the 
LDOS are shown. 



FIG. 3: (Color online) The images of the LDOS with the pure non- 
magnetic potential U = 5.0 eV at different bias voltages. 



stronger than that of the TRS near the impurity. However, 
when co - 120 meV, the LDOS possesses the local minima at 
r = AT mb - 27mi, shown in Fig. 3(c). This means that the 
effect of the TRS on the LDOS is dominant far from the im- 
purity. Therefore, the Friedel oscillations induced by the QPI 
decay faster than those due to the TRS. 

To further elucidate the origin of the dips at the backscat- 
tering wave vectors, we also calculate the LDOS produced by 
a magnetic potential V m = 10.0 eV and U = eV, where the 
TRS is broken. Therefore, the backscatterings are allowed. 



When oj = -120 meV, the LDOS has higher peaks at r = T M2 , 
opposite with a pure nonmagnetic impurity, while there are 
also the dips at r = Tk2- However, p(r, oj) has no obvious 
peaks or dips at qMB and q^#. So the QPI can be neglected 
at the backscattering wave vectors. Therefore, the dips at the 
backscattering wave vectors are indeed produced by the TRS 
in the presence of a nonmagnetic impurity. When the energy 
increases, the peaks and the dips also move forward to the 
magnetic impurity, similar to a nonmagnetic potential. In Fig. 
4(b) and 4(c), the LDOS has the local minima at r = T Ki and 
T M \ , contrary to a nonmagnetic impurity. 
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FIG. 4: (Color online) The images of the LDOS with the pure mag 
netic potential V m = 10.0 eV at different bias voltages. 



In summary, we have proposed a method to detect the TRS 
in topological surface states. Due to this symmetry, the topo- 
logical surface states with opposite momentum and the same 
s are incoherent and contribute zero to the LDOS. Therefore, 
the dips at the backscattering wave vectors in the LDOS in- 
duced by a nonmagnetic impurity can be regarded a signature 
of the TRS. We note that in each modulation direction there is 
the competition between a backscattering wave vector and two 



QPI wave vectors, which leads to robust dip-hump features in 
the LDOS. Because the peaks induced by the QPI decay faster 
than the dips produced by the TRS, it is easy to observe the 
local minima at the backscattering wave vectors at higher bias 
voltages and several periods away from the nonmagnetic im- 
purity by STM experiments. 
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